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Abstract 

We show an equivalence between the approach of Buttiker and the Fermi 
quantum stochastic calculus for mesoscopic systems. To illustrate the method 
we derive the current fluctuations in a two terminal mesoscopic circuit with 
two tunnel barriers containing a single quasi bound state on the well. The 
method enables us to focus on either the incoming/outgoing Fermi fields in 
the leads, or on the irreversible dynamics of the well state itself. The quantum 
stochastic calculus we use is the Fermi analogue of the input /output methods 
of quantum optics. 
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I. INTRODUCTION 



The theory of conductance in mesoscopic electronics was developed some years ago by 
Buttiker, following upon earlier ideas of Landauer H]. The conductance of a mesoscopic sys- 
tem is given in terms of the scattering within and between quantum channels and involves 
the transmission and reflection coefficients as well as the thermal occupation of reservoirs 
feeding or draining those channels. The theory makes direct contact with measured currents 
through averages of quadratic functions of Fermi field operators in the channels. The com- 
putation of the scattering matrices depends on the nature of the systems connected to the 
reservoirs, which could be a simple tunnel barrier or an array of coherently coupled quantum 
dots. In the Buttiker approach once the scattering matrices are calculated we do not need 
to refer to the dynamics of any local systems to which the input and output channels are 
coupled. In many ways this theory is a Fermion analogue of the quantum description of 
optical fields interacting with an optical cavity under the Markov approximation. Such sys- 
tems are described by the input/output theory of Collett and Gardiner 0. The properties 
of the fields outside the cavity are determined by a scattering matrix connecting the input 
and output fields to the cavity and the dynamics inside the cavity. The input/output theory 
for optical fields has been shown to be an example of the quantum stochastic calculus for 
Boson fields 0. 

The Buttiker approach is particularly useful in determining measured properties of the 
mesoscopic system, such as conductance. However recent interest in coherently coupled 
quantum dots for quantum computation [^,|^ has focussed attention on the dynamics of 
localised systems, such as quasibound states on quantum dots, rather than the transport 
through input and output channels. In the input/output theory of quantum optics, the 
dynamics of the local system is described through a Markov master equation, and there is a 
consistency between a description entirely in terms of the input and output modes and the 
irreversible dynamics of the local system to which they couple. In this paper we establish the 
connection between a description in terms of input and output channels and the irreversible 
dynamics of a localised quasibound state on a single quantum dot. The analysis is easily 
extended to more complex local systems. This enables us to make a connection between 
two very successful theories, quantum optics and quantum mesoscopics, for the treatment 
of quantum stochastic processes. We expect that analogies between quantum optics and 
mesoscopic electronics will prove useful as the latter explores the physics of strong coherent 
coupling between local systems (eg quantum dots 0), quantum limited measurements (eg. 
using single electronics [0-0]), and proposals for quantum computation [^0. 

Our treatment will be based on a particularly simple system; a single quantum dot 
coupled to two quantum channels, see figure as this system is the electronic analogue of 
a single Fabry-Perot cavity in quantum optics. The conductance properties of this system 
can easily be obtained by the method of Buttiker M . Recently a similar result was obtained 



using a Markov master equation description of the quasibound state of the dot |TT]]. In 
this paper we give an equivalent description in terms of the quantum stochastic calculus for 
Fermi fields. 
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II. QUANTUM STOCHASTIC CALCULUS FOR FERMIONS 



The free Hamiltonian for a Fermion channel is 

H = hJ2^kaW (1) 

k 

where ak is a Fermi annihilation operator satisfying the anti commutation relations 

aktti + aiUk = (2) 
a^al + a\ak = Ski (3) 



We now define the Fermi field operator 



a(t) = ^^afcC-^^"*-"")* (4) 

k 

which represents free field dynamics with respect to a frame rotating at frequency uq. This 
frequency will later be taken to characterise the energy of the quasi bound state to which the 
Fermi field is coupled. The field operator satisfies the continuum anti commutation relations 

a{t)a\t') + a\t')a{t) = 6{t - t') (5) 

Our objective is to define a quantum stochastic process to accurately characterise the Fermi 
statistics of these fields. To that end we define the integrated operators, 

A{t)= f dt'a{t') (6) 



and the corresponding Ito increment 

dA{t) = A{t + dt) - A{t) (7) 

We now need to specify the state of the free fields. We will take these to be thermal 
equilibrium states of a noninteracting Fermi system at temperature T. It is then easy to 
show, under appropriate assumptions that 

{dA^) = {{dA^f) = (8) 
{dA^dA) = f{uJo)dt (9) 
{dAdA^) = (1 - f{uo))dt (10) 

where the equilibrium occupation probability f{E) is evaluated at the reference energy 
Eq = hujQ. Later this will be the probability that a free field Fermi state, resonant with 
quasi bound state, is occupied. The important point to note here is that these quantities, 
while quadratic in the field increments, are only first order in the time increment. This 
is a quantum analogue of the classical Wiener stochastic process . Stochastic integrals 
of averaged field operators are found using a generalisation of the Ito calculus for classical 
stochastic processes [IT3. In particular we have. 
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/ dt / dt'(dAUt)dA(t'))X(t)Y(t') = / ' dtf(uo)X(t)Y(t] (11) 

JO JO JO 

/ dt / rft'(rfA(t)rfA^(t'))^(t)l^(t') = / ' dt{l - f{uJo))X{t)Y{t) (12) 

JO JO JO 



where X{t),Y{t) may be operator valued functions. We also note that in the Ito calculus 
the chain rule is modified according to 

d{A{t)B{t)) = dA{t)B{t) + A{t)dB{t) + dA{t)dB{t) (13) 

The final term must be included to ensure a correct expansion to linear order in the time 
increment . 

We need now to specify how the free field is coupled to a local electronic degree of 
freedom described by the Fermi annihilation and creation operators c, This degree of 
freedom could for example be a quasi bound state of a single quantum dot or donor atom. 
This interaction will be taken to be linear in the system operators. The usual tunnelling 
interaction is then specified by the Hamiltonian, 

Hi = J2 9kalc + glakC^ (14) 

k 

Such a coupling is the usual way to describe tunnel coupling between a reservoir and a 
localised degree of freedom [|l^] We now anticipate the Markov approximation by assuming 
that around the reference energy, the coupling constants, gk are very slowly varying functions 
of k and replace them by a constant, ^/7. For further discussion see 0. 

Hj = ^ica\t) + c^a{t)) (15) 

The time evolution operator over the time increment dt is given by 

U{dt) = exp{-i^{cdA\t) + cUA{t))} (16) 

This enables us to define the 'output' field stochastic process as 

dAoutit) = U\dt)dA{t)U{dt) (17) 

This equation suggests that we regard dA{t) as the 'input' field stochastic process dA{t) = 
dAin{t). It is then easy to see that 

dAoutit) = dAUt) - Zv^(2/(wo) - l)c{t)dt (18) 

This expression may be written directly in terms of the reservoir field operators as 

aoutit) = ain{t) - i^{2f{uJo) - l)c(t) (19) 

This expression can be used to establish a connection between input and output fields once 
the dynamics of the local system operators is given. In the frequency (energy) domain the 
resulting expression is equivalent to the scattering matrix in the method of Buttiker. 

The dynamics of the local system is specified by a quantum stochastic differential equa- 
tion, for example 
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dc{t) = U\dt)c{t)U{dt) - c{t) (20) 
where the noise operators are 



-lc{t) + i^dH{t) (21) 



dH{t) = dA{t){c{t)c\t) - c\t)c{t)) (22) 

It is at this point that we recognise the difference between Fermi and Bose quantum stochas- 
tic calculus. In the Bose case the noise operator does not depend on system operators, and 
is in fact simply given by dA{t), as the commutator in Eq ^ is unity. In the Fermi case 
however the noise depends on the system operators. This will ensure that the dynamics of 
the local system reflects Fermi statistics, in particular it will ensure that the stochastic dif- 
ferential equation for the number operator nc{t) = c'{t)c{t) takes the correct form. Despite 
the fact that the noise operators depend on the system operators, the average over quadratic 
combinations do not depend on the system. This is because the system operator appearing 
in Eq(p^) has only two eigenvlaues ±1, which when squared is unity. Thus we find 

dH\t)dH{t) = dA\t)dA{t) (23) 

The quantum stochastic differential equation for the system number operator is given by 

dn^it) = c\i)dc{t) + dc\t)c{t) + dc\t)dc{t) (24) 

= -^{n,{t)- f{uj,))dt + dN{t) (25) 

where the noise operator is 

dN{t) = i^{c{tydH{t) - c{t)dH\t)) (26) 

Equation ^ correctly reflects the Fermi statistics of the local system. If the Ito correction 
term, dc\t)dc{t) had been neglected, this would not have occurred. When this term is 
evaluated, the dependence of the noise operators on system variables in Eq. ^ is crucial. 
As a result the average occupation of the local system in the steady state is given correctly 
by 

{ric) = fiooo) (27) 

as would be expected for a Fermi particle. 

We can also obtain the master equation for the local system by 

dp{t) = Tiacen [u{dt)W{t)U^ (dt) - W{t)] (28) 

where W{t) is the state of the total system (localised system and the external fields), and 
Trace refers to a trace over external field variables. Using the noise moments we find, 

^ = -t[H, Pit)] + 7/(a;o)(2ctpc - cc^ p{t) - p{t)cc^) (29) 

+ 7(1 - /(cuo))(2c/9ct - c^cpit) - pit)c^c) 

It is an easy matter to verify that the mean occupation of the dot {nc{t)) obeys the same 
equation that results from taking moments of both sides of Eq. |2^. The master equa- 
tion represents the Schrodinger picture dynamics, while the quantum stochastic differential 
represents the Heisenberg picture. 
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III. QUANTUM STOCHASTIC DYNAMICS OF A SINGLE QUANTUM DOT 



The system we will discuss in this paper is a standard mesoscopic configuration |TB| in 
which ohmic contacts couple to propagating channels to either side of a quantum dot (see 
figure 0). We suppose that there is a single quasi bound state between two tunnel barriers. 
Spin will be ignored. It can easily be included as another state in the dot. We will also 
ignore Coulomb blockade, which for a single bound state simply leads to a shift in the energy 
of the state. We also take the so called 'zero-temperature' limit, and assume that the bound 
state energy is below the effective Fermi energy in the source (L, in figure |l|) and above the 
effective Fermi energy in the drain (R in figure |T]). The ohmic contact at the left of the 
dot is assumed to be a perfect emitter while the ohmic contact at the right of the dot is a 
perfect absorber. This ensures that both reservoirs will remain close to thermal equilibrium 
at all times, provided they are connected to an external EMF. These assumptions will be 
important when we consider the measured quantities in this system. 

We need now to specify how the free fields in the left and right channels are coupled to a 
local electronic degree of freedom described by the Fermi annihilation and creation operators 
c, . This degree of freedom could for example be a quasi bound state of a single quantum 
dot or donor atom. This interaction will be taken to be linear in the system operators. 
The usual tunnelling interaction, anticipating the Markov approximation as in the previous 
section, is then specified by the Hamiltonian, 

Hi = Vll{ca{{t) + c^aL{t)) + v^(cajj(t) + c^aR{t)) (30) 

where 7^, 7j? refers to the tunnelling rate across the left and right barrier respectively, while 
ci'L{'t),ciR{t) are the Fermi fields in the left channel and right channel respectively. The 
corresponding unitary evolution operator for a time increment dt is 

U{dt) = exp [-i^{cdA[ + cUAl) - i^icdA^ji + cUAr)) (31) 

where dA^lt) , dA^lt) are the quantum stochastic processes in the left and right channels 
respectively. The input/output relations are then found to be 

dAL,out{t) = dAL,in{t) - i^c{t)dt (32) 
dAR^out{t) = dAR^in{t) + i^J^c{t)dt (33) 

with the following averages for the noise 

{dAL{t)^dAL{t)) = dt (34) 
{dAR{t)dA^ji{t)) = dt (35) 

All other averages are zero. Alternatively we may write the input-output relations as 

aL,out{t) = daL,in(t) - iy/rLc{t) (36) 

aR,out{t) = ttR^init) + iy/^c(t) (37) 

The quantum stochastic differential equation for the destruction operator in the dot is. 
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dc{t) = -^^l^^^^c{t)dt + i^dHLit) + i^dHRit) (38) 

The quantum stochastic differential equation for the number operator on the dot is 

dnc{t) = 7i(l - nc{t))dt - -iRnc{t)dt + dNiit) + dNuit) . (39) 

If we take moments of both sides of Eq.^ we find, 

djx 

— = 7l(1 -n) - -fun (40) 

The first term corresponds to injection from the source onto the dot. This term is zero if 
the dot is already occupied and n{t) = 1. The second term correspond to emission from the 
dot through the right barrier into the drain. The steady state occupation number on the 
dot is 

noc = — ^ (41) 
The master equation for the dot is found to be 

^ = = Y (2ctpc - cctp - pcct) (42) 



c'cp — pec 



This equation was previously derived by more direct methods in reference [11 



IV. WHAT IS MEASURED ? 

It is at this point we need to make contact with measurable quantities. In the case of 
electron transport, the measurable quantities reduce to current I{t) and voltage V{t). The 
measurement results are a time series of currents and voltages which exhibit both systematic 
and stochastic components. Thus I{t) and voltage V{t) are classical conditional stochastic 
processes, conditioned by the underlying quantum dynamics of the quasi bound state on the 
dot. The reservoirs in the ohmic contacts play a key role in defining the measured quantities 
and ensuring that they are ultimately classical stochastic processes. Transport through the 
dot results in charge fluctuations in either the left or the right channels. These fluctuations 
decay extremely rapidly, ensuring that the channels remain in thermal equilibrium with the 
respective ohmic contacts. For this to be possible charge must be able to flow into and out 
of the channels from an external circuit. We assume that a constant potential difference 
is maintained between the two reservoirs either side of the dot. While the entire system is 
clearly not in thermal equilibrium, we assume that the left and right channels are themselves 
close to thermal equilibrium and can each be specified by a separate chemical potential p^ 
and pr, and these are held constant by a external voltage source, V . 

If a single electron tunnels out of the dot into the right channel between time t and 
t + dt, its energy is momentarily above the Fermi energy. This electron scatters very strongly 
from the electrons in that channel and propagates into the right ohmic contact where it is 
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perfectly absorbed. The net effect is a small current pulse, dliit) in the external circuit. 
This is completely analogous to perfect photodetection: a photon emitted from a cavity will 
be detected with certainty by a detector which is a perfect absorber. Likewise when an 
electron in the right channel tunnels onto the dot, there is a rapid relaxation of this unfilled 
state back to thermal equilibrium as an electron is emitted from the right ohmic contact 
into the depleted channel. This again results in a current pulse in the circuit connected to 
the ohmic contacts. The energy gained when one electron is emitted from the left reservoir 
is, by definition, the chemical potential of that reservoir, fii while the energy lost when 
one electron is absorbed into the right reservoir is hr. The net energy transferred between 
reservoirs is ^l—^r. This energy is supplied by the external EMF, V and thus fiL—fJ'R = eV- 
It should not be supposed that the electron injected from the left contact and emitted into 
the right contact have the same energy as the energy of an electron on the dot. In fact any 
electron energy at all will suffice to restore thermal equilibrium in the left and right leads. 
If an electron is emitted into the left channel between times t and t + dt, the (unnormalised) 
sate of that channel is PL,o{t + dt) = a\ Qp{t)aL^odt. The probability of this event occurring 

is simply the normalisation constant and is Pe{t) = ii^icLL oPi^)^L,o)dt = {aL,o{t)(iL oi^)) ■ 
That is to say the probability of emission of electrons into the left channel is determined 
by the anti normally ordered number flux operator in the left channel. A similar argument 
indicates that the probability of absorption of an electron in the right ohmic contact is 
given by the mean of the normally ordered number flux operator in the right most channel, 
Pait) = {0'Roi^)^R,o{t)) ■ This is precisely analogous to perfect photodetection from an 
optical source Both the emission into the left channel and absorption from the right 
channel are conditional point processes, conditioned on the quantum state of the quasi bound 
state on the dot. 

On average of course the same current flows in both reservoirs, however as the current is 
stochastic it is made up of contributions from pulses in each lead, which do not necessarily 
occur at the same time. Indeed they are necessarily separated in time by a degree depending 
on the life time of the quasi bound state in the dot. In mesoscopic devices however current 
in measured locally in each lead, thus we can consider either the current in the left lead, 
lL{t) or the current in the right lead and correlations between them. The current that 
flows in the right lead is simply given by the probability per unit time that an electron in 
that channel is absorbed by the perfect absorber that is the right Ohmic contact. Thus 



The average current that flows in the left lead is given by the average probability per unit 
time that an electron is emitted by the perfect emitter that is the left Ohmic contact. 



Note that the average on the left hand side is an average of a classical stochastic process, 
while the average on the right is of a quantum stochastic process. We may now substitute 
the relationship between the output flelds, the input flelds and the operator for the local 
state, Eqs(37,38 ). The average currents are then found to be 



E{Ijt{t)) = e(a}j,oWai?,o(t)) 



(43) 



E{h{t)) = e{aL,oityL,o{t)) 



(44) 



Eihit)) = e7L(l - Mm 



(45) 
(46) 
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In the stationary state both currents are equal and given by 



(47) 
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where 7 = 7l + Tr- 

The stationary two time correlation matrix is given by 

GaAr) = E{Ut + r), Ip{t))t^oo (48) 

and E{X,Y) = E{XY) - E{X)E{Y). The quantity E{I^{t + T)Ip{t) is determined by the 
appropriately ordered two time correlation function for the quantum fields in the channels. 
As both emission and absorption are point processes we find, 

Eihit + r)/i(t) = e^aL,o{t)ai^oit))^i^) + e^aL,oit)aL,o{t + T)al o(t + r)4 oW)r>o (49) 

E{lR{t + T)lR{t)) = e^aR^o{tyaR^o{t))6{r) + e^{a^R^oit)aR,oit + ^)aR,o{t + r)aR,o{t))r>o (50) 

E{lR{t + r)/i(t)) = (aJj,o(^ + + ^)aL,o(t)4,oW) (51) 

E{h{t + r)/^(t)) = {aL,o{t + T)ai,o(t + r)aJj,o WaR,o(t)) (52) 

Using equations 37 and 38, the steady state {t 00) field correlation functions may be 
expressed solely in terms of the correlation functions for the quasibound state as 

{aL,o{t)aL,o{t + T)a{ o{t + 'r)al,o(0)r>o = 7LTr (^cc^e'^^cVooc) (53) 

{(^R,oit)(^R,oit + 'r)aR^o{t + T)aR,o{t))T>o = 7^Tr {^c^ ce^^'' cp^oC^^ (54) 

(«ko(^ + + '^)«i^,o(^)«L,o(0) = 7i?7LTr (^c^ce^^cVooc) (55) 

{aL,o{t + T)aL.o(^ + T)(^R.oit)"-R,o{t)) = iLlR^'c (cch^'^ cpooC^) (56) 



where poo is the steady state solution to the master equation, Eq ^ for the quasi bound 
state on the dot. 

Upon solving the master equation we may evaluate each correlation function to give; 

E{h{t + r), hit) = E{lR{t + T)jR{t) = e2^5(r) - e^^e"^^ (57) 

E{lR{t + r), /z.(t) = e^^e--^- (58) 

i?(/i(t + r), /^(t) = e^^-^e-^- (59) 

where 7 = 7l + 7_r. The power spectrum of the noise in either the left or the right lead is 
then given by the Fourier transform of the current correlation function in either lead 

Sr{uj) = Sl{uj) = eioo I 1 - ] (60) 

Note that at zero frequency and symmetric rates (7/? = 7^) the current noise is suppressed 
by a factor of 0.5 over shot noise. This is the same as that obtained by Buttiker [|10|, and 



is equivalent to that obtained in reference [11 1. Note however that in that paper it was 



assumed, following earlier work |T^, that the measured current is a superposition of the two 
Poisson processes of emission and absorption. This however is not the case for mesoscopic 
measurements. As has been stressed by Buttiker, the current is measured in one lead at a 
time and thus the correct expression for the current noise is that given above. 
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V. DISCUSSION AND CONCLUSION 



In this paper we have shown an equivalence between the approach of Buttiker and the 
approach of quantum stochastic calculus to the current through mesoscopic systems. To 
illustrate the equivalence we have discussed the current fluctuations in two terminal meso- 
scopic circuit with two tunnel barriers containing a single quasi bound state on the well. 
The method enables us to focus on either the incoming and outgoing Fermi fields in the 
leads, or on the irreversible dynamics of the well state itself. We have of course made the 
Markov approximation in order to obtain the quantum Langevin equations for the system. 
The Markov assumption is equivalent to the assumption of the Breit-Wigner (Lorentzian) 
assumption for the transmission coefficients though a double barrier structure, and is dis- 
cussed in some detail in reference |TT|. 



We believe there are two advantages in our approach. Firstly it is useful to be able to 
refer to the quantum irreversible dynamics of the quasi bound states on local systems defined 
by barriers, as well as the input and output Fermi fields. This is particularly important in 
coherently coupled quasibound states. This is essential for recent condensed matter schemes 
for quantum computation, where the focus is not so much on the properties of the external, 
classical currents but rather on the dynamics of the local systems themselves. Secondly, our 
method parallels a similar approach to strongly coupled field modes in quantum optics, thus 
suggesting useful directions for future work. We mention one such direction. In quantum 
optics the method of quantum trajectories |TB[ enables a description to be given of quantum 
limited measurement ||T^, quantum control (feedback and feed forward) pO[, and cascaded 
local systems (irreversible but directional coupling) |2l| . As mesoscopic technology advances 
these topics will become increasingly important. The presentation in this paper shows how 
the quantum stochastic methods of quantum optics may be taken over to Fermi systems. 
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FIGURES 

FIG. 1. Schematic representation of tunneling through a single quantum dot. 
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